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Abstract 

We prove a new theorem about the relationship between optical field Wigner function's 
Radon transform and optical Fresnel transform of the field, i.e., when an input field ip (a/) 
propagates through an optical [D (-B) (— C) A] system, the energy density of the output field 
is equal to the Radon transform of the Wigner function of the input field, where the Radon 
transform parameters are D,B. We prove this theorem in both spatial-domain and frequency- 
domain. 

In optical communication theory every signal or image can be uniquely and indirectly described 
by a Wigner distribution function (WDF) [U [21 [3] . The WDF (or named Wigner transform) of an 
optical signal field ip (x 1 ) is defined as 
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Wtp(v',x') involves both spatial distribution information and space-frequency distribution informa- 
tion of the signal, v is named space frequency. W^{v' ', x') is said to be bilinear in the signal because 
the signal enters twice in its definition. The WDF undergoes certain variations if something happens 
to the signal. For examples, passage through a lens corresponds to a vertical shearing of the WDF, 
propagation in free space means a horizontal shearing of the WDF [3] . However, the WDF preserves 
space and space frequency marginal properties of any signal, 
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dv'W^v' ,x') = \ip{x') | 2 , (2) 
dx'W^v' ,x') = \4>(v')\ 2 , (3) 



where ip {v) = f_™ ~^§^' t P ( x ) e%Xl ' ■ ^ one wan ts to reconstruct the Wigner function by using various 
probability distribution, obviously the position density \ip (x 1 ) | 2 and the space-frequency density 
\tj) (V) | 2 were not enough, so the Radon transform [51 [3] of the Wigner function is introduced [7], 
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R (x) = {{ dx'dv'8 (x - Dx' + Bu') W^{v' , x'), (4) 



R (x) is also a probability distribution along an infinitely thin phase space strip denoted by the 
real parameters D, B. The inverse relation of ([4]) is the foundation of optical tomographic imaging 
techniques (the techniques derive two-dimensional data from a three-dimensional object to obtain a 
slice image of the internal structure and thus have the ability to peer inside the object noninvasively.) 



On the other hand, an optical system can be analyzed by using either diffraction theory or ray 
optics, ft is worth setting up a connection between ray transfer matrix and diffraction theory, i.e., the 
diffraction integration written in terms of ray transfer matrix — named Collins' diffraction integral 
formula [HI EE HO] which describes the propagation of a general beam ip {x') through an {ABCD) 
optical paraxial system, 
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where (p (x) is the output field, AD — BC = 1, and K. (x, x') is the integral kernel 
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If the energy in the initial beam is normalized, J^°° \ip (x')\ 2 dx' = f, then the output beam is 

normalized too, J^°° \(p (x)\ 2 dx = 1. Clearly, if the [ABCD] system is changed to [D (~B) (— C) A] 
system, then Eq. (JSJ) should read 
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where the integral kernel K. (x, x') is 
JC (x, x') 



IC {x, x') ip {x') dx' , 
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In this Letter, we shall reveal the following theorem: 

When an input field ip (x') propagates through an optical [D (~B) (— C) A] system, the 
energy density of the output field </> (x) is equal to the Radon transform of the Wigner 
function of the input field, where the Radon transform parameters are D, B. The proof 
is demonstrated as follows. 

With the use of Dirac 5-function, we can re-express the Wigner function of input field ip (x') in 
Eq. (TTD as 
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Substituting ([9]) into Q we rewrite the Radon transform of W^iy 1 , x r ) as 
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On the other hand, when the beam ip [x') propagates through the [D (-B) (— C) A] optical system, 



according to the Fresnel integration (O-©, we have 

\<f>{x)\ 2 = / dx'iC{x,x')^{x') / dx"ip*(x")iC*(x,x") 
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which is the same as R(x) in |10]). So combining g]), dTO])- (HU) , and © we reach the conclusion 
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(12) 

where DA — BC = 1. The physical meaning of Eq. (fT2"]) is: when an input held propagates through 
an optical [D (— B) (— C) A] system, the energy density of the output field is equal to the Radon 
transform of the Wigner function of the input field. So far as our knowledge is concerned, this 
conclusion seems new. 

Eq. ([8]) is the relationship between the input amplitude and output one in spatial-domain. Now 
we turn the above discussion to the case of space-frequency domain. 

For a [D (-B) (— C) A] optical system, the Collins' diffraction integral formula in space-frequency 
(angle spectrum pj]) domain is p~2] 
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where the kernel is 
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(13) 
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In space-frequency domain the Radon transform is along an infinitely thin phase space strip denoted 
by the real parameters A, C, 



On the other hand, the Wigner function expressed in terms of the space-frequency field is 



R{v)= // dx'dp'5(u-Av' + Cx')Wr(u',x') 



(16) 



rather than the parameters D, B of the Radon transform (j4]) in spatial-domain. We want to examine 
if the above conclusion still holds in the space-frequency domain. For this propose, we rewrite (|15p 
as 
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Then we substitute (fT7|) into (fTBj) . the result is 
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dv"dv'"r (v") $ (v'") exp <j [A (v" 2 - v'" 2 ) - 2v iy" - v'")\ \ . (18) 



On the other hand, from (11311 we calculate 
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Comparing Eq. |19]) with Eq. (18]) and using ([I3|)-([l4]) we see 
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this is the theorem expressed in space-frequency domain. 

We now take an example to confirm the theorem. When the input field is described by a Gaussian 
function (Gaussian chirplet) 

(21) 
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which depicts a Gaussian windowed linear chirp signal [13] • Using Eq.fT]) we obtain its Wigner 
function 
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which shows that the energy of Gaussian chirplet is concentrated at v' — 0x' . According to Eq. (J4|, 
the Radon transform of (|22|) with the parameters D, B is 
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On the other hand, according to and ([8|), the Fresnel transform of input field W$ (v\ x') through 
a [D (-B) (-C) A] optical system is given by 
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Comparing (12"4")) with (|2"3")) and noticing 
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we see 

|0o (x) | 2 =#o(z), (26) 

as expected. 

In summary, we have derived a new theorem governing the connection between optical field's 
Fresnel transform and its Wigner function's Radon transformation, since both the Wigner function 
and the Fresnel transform are widely used in optical propagation, we hope this theorem would 
have new applications in the analysis of optical communication and optical tomography. For the 
application of Wigner function in deriving the quantum-mechanical photocount formula, we refer to 
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